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1.  EXECUTIVE  SUMMARY 


Statement  of  the  Problem  Studied 

This  seedling  study  was  meant  to  evaluate  the  potential  of  magneto-electric  materials 
to  impact  size,  weight  and  performance  of  tactical  radio  and  radar  systems.  ME  materials 
are  a  subclass  of  multifunctional  materials  that  are  distinguished  by  simultaneously 
having  two  or  more  of  the  following  characteristics;  ferromagnetism,  ferroelectricity  or 
ferroelasticity.  As  a  consequence  of  this  combined  functionality  such  materials  can 
exhibit: 

•  Wide  dynamic  frequency  response  as  a  consequence  of  spontaneous  polarization 
and  spontaneous  magnetization 

•  Frequency  dependent  variable  penneability  and  frequency  dependent  variable 
permittivity  in  a  single  device 

•  Magnetic  properties  controlled  by  electric  field. 

The  goals  of  the  seedling  were  to: 

•  Compile  a  database  of  material  parameters  specific  to  ME  materials. 

•  Perform  device  simulations  for  a  tunable  filter  architecture  using  materials 
parameters  from  the  ME  database  for  input. 

•  Benchmark  tunable  filter  ME  devices  against  state-of-the  art  devices  presently  in 
use  for  tactical  radio  and  radar. 

•  Solicit  expert  input  on  risk/benefit  analysis  of  tactical  radio  and  radar  systems 
based  on  ME  devices. 

Summary  of  Results 

A  database  of  properties  for  known  single  phase  and  composite  ME  materials  was 
compiled  which  clearly  showed  that  only  composite  ME  materials  were  suitable  at 
present  for  technological  applications.  A  parametric  study  was  performed  using  a 
heuristic  model  to  show  qualitatively  the  effect  of  various  materials  parameters  on  device 
tunability. 

The  benchmark  ME  device  was  a  composite  of  a  piezoelectric  and  magnetostrictive 
ferrite  which  displayed  tunablilty  of  the  ferromagnetic  resonance  by  means  of  strain 
coupling  (a  second  order  ME  effect).  A  significant  impediment  to  achieving  our  goals  of 
modeling  this  device  was  the  discovery  that  existing  commercial  codes  which  claimed  to 
treat  ferrite  materials  could  not  accommodate  the  nonuniform  magnetic  fields  associated 
with  common  device  geometries  and  could  also  not  incorporate  magnetostriction  in  any 
meaningful  way.  Development  of  alternative  code  was  beyond  the  scope  of  this  seedling. 
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2.  INTRODUCTION 


Modern  military  tactical  radio  requires  the  ability  to  receive  and  transmit  over  an 
extremely  large  band  of  frequencies,  from  the  low  MHz  well  into  the  GHz  range.  This 
broadband  environment  differentiates  military  systems  from  the  narrow  bandwidth 
requirements  of  commercial  technologies.  The  need  for  simultaneous  regard  over  a  very 
large  bandwidth  places  very  great  demand  on  the  linearity  of  the  active  and  passive  signal 
processing  in  military  communication  systems.  Improvements  in  passive  filter 
technology  are  expected  to  relax  the  demands  on  the  active  electronics  in  these  systems, 
but  at  the  expense  of  an  impractical  number  of  fixed  frequency  filters,  or  through  the  use 
of  a  more  reasonable  number  of  tunable  filters.  Varactor  tuned  filters  can  provide  the 
required  tuning  range,  but  varactors  are  inherently  non-linear  due  to  the  solid  state 
physics  behind  the  tuning  mechanism.  A  key  technology  for  flexible,  low-power,  and 
compact  tactical  radios  is  improved  tuning  elements  for  passive  filters.  Military  radar  has 
similar  needs. 

2.1  Issues  Potentially  Addressed  by  ME  Materials  and  Devices 

Like  most  applications,  engineers  of  microwave  and  millimeter  wave  devices  and 
subsystems  trade  cost,  design  and  performance.  When  evaluating  different  technologies 
they  look  for: 

•  Component  count  vs.  frequency, 

•  Tunability,  cost  vs.  hysteretic  loss, 

•  Cost  vs.  frequency  tunability, 

•  Frequency  tunability  vs.  hysteretic  loss, 

•  Resistivity  vs.  quality  factor, 

•  Degrees  phase  shift  /  db  loss  vs.  tunability, 

•  Insertion  loss. 

Specifications  that  are  important  for  nonreciprocal  devices,  such  as  circulators  and 
isolators,  are  saturation  magnetization,  easy-axis  or  easy-plane  anisotropy,  anisotropy 
magnitude,  magnetoelastic  effects,  dielectric  constant,  resistivity,  power  handling 
capability  (limited  by  spin  wave  instabilities).  These  material  properties  can  be 
engineered  to  meet  the  designer’s  performance  criteria  above  for  many  RF  applications. 
Y et  increasing  communications  demands  and  needs  have  created  a  technological  pull  for 
materials  and  devices  which  require  higher  performance.  Some  of  the  “game  changing” 
device  characteristics  brought  up  by  the  RF  engineering  community  are: 

•  Circulators 

o  Control  impedance  for  improved  uniformity; 

•  Phase  shifters 

o  True  time  delay 

o  Independent  control  of  phase  and  amplitude; 

•  Power  dividers  /  combiners 

o  Impedance  transfonning 

o  Impedance  maintenance; 

•  Low  loss,  high  s,  wide  band  spiral  radiators  to  protect  LNAs. 
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The  hope  is  that  multiferroic  and  magnetoelectric  materials  and  devices  designed 
from  them  may  help  meet  some  of  challenges  RF  engineers  face.  The  results  of  this 
seedling  show  that  the  understanding  of  multiferroic  materials,  their  capabilities  / 
limitations  and  maturity  level  is  warrants  further  investment. 

3.  TASK  1:  DATABASE  OF  MAGNETO-ELECTRIC  MATERIALS  PROPERTIES 

The  benefits  of  incorporating  magneto-electric  (ME)  materials  in  electronic  devices 
has  been  reinvigorated  by  the  promising  developments  of  monolithic  materials.  Recently, 
intensive  research  has  been  conducted  to  optimize  the  magneto-electric  effect  in  a 
number  of  different  materials  systems.  Beyond  the  initial  discovery  of  inherent  magneto¬ 
electric  effect  in  C^Cfi,  other  monolithic  materials  have  been  discovered  that  have  exhibit 
a  larger  magneto-electric  effect.  Fiebig1  gives  a  nice  review  of  other  ME  materials  that 
have  been  found  and  their  ME  effects:  TfiCfi,  ferrites  such  as  (GaFeCfi,  BiFeCfi  and  YIG), 
boracites  (such  as  Ni3B70i3I),  BaMnF4,  TbP04,  mixed  perovskites  (such  as 
PbFe0.5Nb0.jCfi),  simple  perovskites  (such  as  FIoMnCfi)  and  composites  (such  as  PZT- 
Terfenol-D).  Other  possible  candidate  materials  are  thought  to  include  colossal 
magnetoresistive  oxides,  pyrochlores,  and  other  ferromagnetic  oxides.  Suitability  of  these 
materials  for  application  has  several  criteria. 

Primary  screening  of  ME  materials’  readiness  for  application  is  determined  by 
ferroelectric  and  magnetic  ordering  temperatures.  Both  Curie  and  Neel  temperatures  must 
be  above  the  range  of  interest,  in  our  case  200°C  or  greater.  Secondarily,  the  ME 
coefficient  must  be  a  sufficiently  large  to  be  of  technological  interest.  Thermodynamic 
limits  of  the  linear  ME  tensor  in  monolithic  materials  are  determined  by  electric,  y? ,  and 
magnetic  susceptibilities,  f  \  according  to  a2 3i}  <  XuZ'jj  •>  where  ay  is  the  induction  of 
polarization  by  a  magnetic  field  (direct  effect)  or  magnetization  by  an  electric  field 
(converse  effect).  In  composites,  an  effective  ME  effect  occurs  as  a  product  property, 
which  refers  to  novel  effects  stemming  from  the  interaction  between  the  composite 
constituent  materials.  Harshe  et  al.  “  developed  expressions  for  an  effective  ME 
coefficient,  ay  for  several  conditions  based  upon  the  constituent  materials  properties  of  a 
magnetostrictive  phase  and  piezoelectric  or  electrostrictive  phase.  These  expressions  have 
been  revisited  and  refined  in  recent  years  by  Nan  et  al? ,  Bichurin  et  al4 5,  Srinivasan  et 
al.  5  .  The  expression  can  be  accomplished  by  combining  magnetostrictive  and 
piezoelectric  materials  for  example: 

, ....  magnetic  mechanical 

ME  effect  =  - - - x - . 

mechanical  electrical 

Ideal  ME  material  characteristics  include  low  loss,  phase  stability  over  operating 
temperature,  mechanical,  magnetic  and  electrical  ranges  of  interest,  easily  magnetized  (//, 

1  Fiebig,  J.  Phys.  D:  Appl.  Phys.  38  (2005)  R123-R152 

2  G.  Harshe,  J.P.  Dougherty  and  R.E.  Newnham,  Internatl.  J.  Appl.  Electromagn.  Matls.,  v4,  pl45-159, 
(1993) 

3  C.W.  Nan,  G.  Liu,  Y.  Lin,  PRL,  197203,  2005 

4  M.I.  Bichurin,  LA.  Kornev,  V.M.  Petrov,  A.S.  Tatarenko,  Y.V.  Killiba  and  G.  Srinivasan,  PRB,  64 
094409,  (2001) 

5  G.  Srinivasan,  M.I.  Bichurin  and  V.M.  Petrov,  PRB  054402,  2003 
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>  100)  and  polarized  (sr  >  100)  with  low  anisotropy  and  low  coercivity,  high  saturation 
polarization  ( Ps  >  2  Op  C/cm  ),  minimal  defects  for  hysteretic  scattering  and  oxygen 
diffusion  and  being  amenable  to  a  scalable  manufacturing  process.  To  date,  most 
literature  focuses  on  the  direct  ME  effect,  while  of  great  interest  for  radio  and  radar  is  the 
converse  effect,  which  has  been  illustrated  at  resonance  by  Shastry  et  al 6  and 
quasistatic  ally  by  Ramesh7 8.  Of  importance,  Ramesh  reported  coupling  coefficient  of 
1x10'"  G  cm/V.  Enhanced  properties  can  be  achieved  by  exploiting  epitaxial  strains 
between  thin  films  and  device  engineering  for  near  resonance  operation  to  maximize  the 
ME  effect. 

The  aim  of  this  proposed  program  was  to  relate  the  unique  characteristics  of  magneto¬ 
electric  materials  to  device,  and  ultimately  system  performance.  The  specific  application 
determines  the  needed  characteristics  of  the  device  materials.  Table  1  presents 
characteristics  of  candidate  ME  materials. 


Table  1.  Performance  and  Characteristics  of  ME  Materials. 


Material 

Current  characteristics 

Reference 

Single  material 

Temperature 

aE  =  4.1  ns/m 

Boracites:  M3B7013X 

M  =  Cr,  Mn,  Fe,  Co,  Cu,  Ni.  X=  Cl,  Br,  1 

er'e0  aE  =0.12  V/cm  Oe,  FE:  Tc~1 00K 

AF  or  weak  FM:  High  coercivity,  saturation 

Ascher  et  al  1966 

Ba/WF4  /  LiMP04 

M=  Mg,  Mn,  Fe,  Co,  Ni,  Zn,  Tb 

aE  =  30-36  ns/m,  FE:  TC~100K 

AF  or  weak  FM:  H  =  20T,  TN  ~25  K 

Scott  1979 

Rivera  1994, 

Kornev  et  al  2000 

Perovskites:  AB'^B"^ 

A  =  Ba,  Sr,  Pb,  B'=  Fe,  Co,  Ni,  B"  =  Cr,  Mo,  W 

Low  transition  temps,  AF  or  weak  FM 

High  fields  necessary 

Manganites:  RMnOs 

R  =  Bi,  Y,  Nd,  Tb,  Dy,  Ho,  Er,  In,  Lu,  Sc,  Tm 
(CMR:  PrCa,  LaSr) 

AF:  TW~70-130K,  High  H,  FE:  7C~1000K 

Rich  mag.  behavior  b/c  frustrated  system 

Kimura  LANL 

Pyrochlores:  A20-B206 

A  =  BiZn,  CaTi,  BiFe,  B  =  ZnNb,  NbTi,  FeNb 

£r'£0  aE=  3  V/cm  Oe,  RT  operation,  high  e,  Magnetic 
variant  of  microwave  dielectrics 

Vanderah  NIST, 
Ubic  et  al. 

Ferrites:  RFe03  (YIG) 

R  =  Ba,  Bi,  Ni,  Co,  Zn,  Mn, 

aE  =  30  ns/m,  Low  resistivity,  high  fields 

Multi-material  fabrication 

Not  scalable,  strain  coupling  dependent,  field  nonuniformities 

Piezoelectric  /  Ferrite 

PZT,  BaTi03  /  Co,  Ni,  Li,  NiZn  Ferrites 

Interdiffusion,  bonding 
aM  acheived! 

Srinivasan  et  al 
2004 

Ryu  et  al  2002 

PZT-PVDF  /  Terfenol-D 

PMN-PT  /  Terfenol-D 

£r'£0  aE  =  4.68,  10.3  V/cm  Oe,  Low  resistivity,  bonding 

Bichurin  et  al  2003 
Ryu  et  al  2002 

4.  TASK  2:  SIMULATIONS 

4.1  Parametric  Analysis  of  Impact  of  Strain  on  Tunability  of  FMR 

o 

Petrakovskii  et  al  have  shown  that  the  frequency  shift  of  the  ferromagnetic 
resonance  for  a  spherical  anisotropic  ferromagnet  subjected  to  a  homogeneous  uniaxial 

6  Shastry,  G.  Srinivasan,  M.  I.  Bichurin,  V.  M.  Petrov,  and  A.  S.  Tatarenko,  Physical  Review  B  70,  064416 
(2004) 

7  R.  Ramesh,  DARPA  DSO  Multiferroics  Workshop,  July  2005 

8  G.  A.  Petrakovskii  and  E.  M.  Smokotin,  Soviet  Physics  JETP  28,  1101  (1969) 
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compressive  stress  in  the  [110]  axis  and  magnetized  by  a  field  H  in  the  (110)  plane  is 
given  by  the  relation; 

1 S  1 S 

+Tn-(*m  -  4„X1  +  cos!(0+  SJ)-2Am)  (1) 

2  Ms 

J _ 4/Tq _ [  Ha  [% sin2  2{0+  8)  — y2  sin2  (61  +  ^)]|  1 

1  Ha  [%  sin2  2(9+  S)  —  y2  sin2  (6*+  fi)]  2H0  j 

Where,  Ho  =  co/\y\  is  the  applied  magnetic  field  in  Oersted,  Ha  =  2K/M  is  the  anisotropy 
field  in  Oersted,  M  =  the  saturation  magnetization  in  Gauss,  6  =  the  angle  between  the 
magnetizing  field  H  and  the  [100]  direction  in  radians,  ^ioo,  =  the  magnetostrictive 
coefficients,  Ki  =  uniaxial  anisotropy  constant,  a  =  the  applied  stress  in  MPa,  and  5  = 
correction  of  the  field  direction  due  to  presence  of  uniaxial  anisotropy  in  radians  (and/or 
elastic  stresses9)  given  as: 

Ha  sin2fl[l  -  %sin2  61] 

2  H0  +  Ha[ 2-  sin2  0-3  sin2  26»]  ’ 

By  inspection,  one  can  see  that  the  magnetic  field  shift  is  linear  with  stress  and 
magnetostriction  coefficients,  while  inversely  proportional  to  saturation  magnetization. 
The  complex  relationship  of  field  angle  and  anisotropy,  given  by  S,  gives  the  correction 
to  theta  due  to  the  magnitude  of  the  magnetocrystalline  anisotropy  as  shown  Figure  1  for 
three  ferrites  of  technological  interest,  YIG,  CoFe2C>4  and  BaM. 


Figure  1.  (a)  Theta  correction  due  to  anisotropy,  (delta)  vs.  theta,  (the  cylindrical  angle  H0  makes 
with  the  |100)-direetion)  for  YIG  CoFe204  and  BaM.  For  most  materials,  this  is  a  small  correction  on 
theta,  (b)  Anisotropy  correction  term  (braces  in  Error!  Reference  source  not  found.)  vs.  theta,  the 
cylindrical  angle  H0  makes  with  the  1 100] -direction  for  YIG  CoFe204  and  BaM.  For  most  materials, 
this  is  a  small  correction  to  5H. 


9  G.  A.  Petrakovskii  and  Yu.  N.  Kotyukov,  Soviet  Physics-Solid  State,  v7,  n8,  (1966). 
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Equation  1  has  been  used  to  model  material  tunability  by  varying  applied  stress  and 
anisotropy-to-magnetostriction  ratio.  Values,  in  Table  2,  for  YIG  and  BaM  are  used 
because  of  technological  importance  and  CoFe2C>4  because  of  its  large  magnetostriction. 


Table  2.  Properties  of  Selected  Ferrites  used  as  input  for  Equation  1. 


Property 

YIG 

CoFe204 

BaM 

MJG1 

1750 

4400 

4420 

H  [Oel 

3400 

12000 

7500 

K  [erg/cc] 

104 

2.6  106 

3.2  106 

^100,^111  [ppm] 

-1.4, -1.6 

-670,  -120 

-13,0 

Figure  2.  Ferrite  tunability  as  a  function  of  theta,  the  cylindrical  angle  FI0  makes  with  the  [100]- 
direction,  at  stress  -20  MPa.  Thus,  maximum  tunability  occurs  when  the  field  is  aligned  with  the 
|100]-direction. 


Figure  3.  Ferrite  tunability  as  a  function  of  applied  stress.  Large  magnetostriction  from  CoFe204 
demonstrates  the  large  effect  that  stress  has  on  tunability.  YIG  has  better  tunability  than  BaM 
because  of  the  much  lower  Ki  values. 


6 


Figure  4.  Ferrite  tunability  as  a  function  of  Kj.  If  K]  were  to  change  from  its  nominal  value  (Table  2), 
then  tunability  increases  as  anisotropy  decreases. 


4.2  Stress  Induced  FMR  Tunability  of  Microwave  Devices 

Theoretical  and  experimental  work  that  investigates  tunability  of  FMR  due  to  stress- 
induced  effects  is  sparse.  Dionne  and  Oates  10  use  a  mechanical  fixture  to  apply 
compression  into  the  plane  of  the  Ni-Al  spinel  ferrite  slab.  They  achieved  tunability  of 
5%  at  10.2  GHz  due  to  high  magnetostriction  with  a  moderate  magnetocrystalline 
anisotropy.  For  YIG  under  compression,  they  achieved  ~50  MHz  frequency  shift  at 
2.9  GHz  resulting  in  1.7%  tunability,  which  they  attribute  to  the  correspondingly  lower 
magnetostriction- to-anisotropy  ratio.  For  a  piezoelectric-controlled  ferrite  device  they 
believe: 

"To  reduce  this  concept  to  practice,  the  transmission  of  stress  from  piezoelectric  to 
ferrite  in  a  planar  configuration  needed  to  control  a  stripe  domain  pattern  must  be 
studied,  to  determine  the  most  efficient  and  practical  method.  Because  the  ferrite  must 
have  a  high  magnetostriction-to-anisotropy  ratio,  the  chemical  composition  of  the 
ferrite  must  be  tailored  to  optimize  that  property  without  degrading  the  permeability 
and  propagation  constants. " 

In  their  patent 1 1  ,  they  suggest  ferrites  that  have  “strong  inverse  magnetostrictive 
properties  sufficient  to  align  magnetic  domains  and  conventional  piezoelectric  materials 
which  impart  mechanical  stress  to  the  ferrite  are  the  most  promising  to  achieve  a  high 
magnetostriction-to-anisotropy  ratio.”  These  compositions  follow: 


10  G.F.  Dionne  and  D.E.  Oates,  Mater.  Res.  Soc.  Symp.  Proc.  v833,  (2005) 

11  G.F.  Dionne  and  D.E.  Oates,  Tunable  microwave  magnetic  devices,  US  patent  6,919,783,  (2005) 


7 


•  YIG  garnet  families 

o  Y  3F  e5-x-yAlxIny0 1 2 
o  Y  3F  e5-x-yGaxIny0 12 
o  Y3F  e5-x-yAlxS  cy0 1 2 
o  Y  3F  e5-x-yGaxScy0 12 

•  Calcium-vanadium  garnet  families 

o  Y3_2XCa2xFe5-x-yVxInyOi2 
o  Y3.2XCa2XFe5-x.yVxScyOi2 
o  Y3.2X.yCa2X+yFe5.x.yYxZryOi2 

•  Lithium,  nickel,  manganese,  and  magnesium  spinel  ferrite  families 

o  Li05+t/2Fe2.5-3t/2TitO4 
o  Lio.5-z/2ZnzFe2.5-z/204 

o  Ni  1  _zZnzF  e2-xAlx04 
o  Mni_zZnzFe2-xAlx04 
o  MgiFe2-xAlx04 

12 

Srinivasan  et  al.  accomplished  a  resonance  field  shift  of  40  Oe  at  E-fields  of 
8  kV/cm  in  a  bilayer  of  4.9pm  (111)  YIG  and  (001)  PMN-PT  with  static  E-  and  H-fields 
parallel  to  <11 1>.  This  is  equivalent  to  112  MHz  shift,  or  1.2%  tunability  at  9.3  GHz. 
The  insertion  losses  from  single-resonator  filters  can  be  theoretically  estimated  to  be  - 
2.5  dB  and  have  a  tunability  of  4.5%  at  100  kV/cm,  though  practically  it  has  been  shown 
to  be  much  greater  than  that,  greater  than  -20  db. 

Theoretical  and  experimental  work  that  investigates  ferromagnetic  resonance  and 
magneto-static  waves  under  non-uniform  conditions  are  few.  Dikshteyn  and  Mal’tsev14 
bent  a  GGG  substrate  with  a  thin  YIG  film  well  below  its  elastic  limit  (0.02%  strain,  non- 
uniform  stress  of  from  56  MPa  at  the  root  to  average  28  MPa)  and  achieved  29  MHz 
frequency  shift  at  a  center  frequency  of  3.5  GHz,  yielding  a  tunability  of  0.8%.  This  is  in 
close  agreement  with  the  results  from  Error!  Reference  source  not  found,  in  Figure  3. 

4.3  HFSS  Simulations  of  Tunable  Filter 

HFSS  (High  Frequency  Structure  Simulator)  is  a  3D  full  wave  finite  element  field 
simulator  code  available  from  Ansoft  Corporation.  The  current  version  of  HFSS  includes 
gyro  tropic  ferrite  materials.  However,  during  the  course  of  this  seedling  it  became  clear 
that  HFSS  has  two  significant  drawbacks  which  were  not  sunnountable  within  the  scope 
of  this  seedling.  First,  HFSS  is  only  capable  of  simulating  problems  with  homogeneous 
magnetic  fields  and  is  only  strictly  applicable  to  an  ellipsoidal  ferrite.  Second,  their 
appears  to  be  no  intuitive  way  to  incorporate  magnetostriction  into  the  model  using  the 
available  inputs  to  the  code. 


12  G.  Srinivasan,  A.S.  Tatarenko  and  M.I.  Bichurin,  Electronic  Letters,  v41,  nlO,  (2005) 

13  P  Carter,  IEEE  Trans.  Microw.  Theory  Tech.,  vl3,  n3,  (1965) 

14  I. Ye.  Dikshteyn  and  O.A.  Mal’tsev,  J.  of  Comm.  Tech,  and  Electronics,  v38,  n5,  (1993) 


4.4  Calculations  of  FMR  as  a  Function  of  Internal  Field 


The  prototype  device  was  a  PZT/YIG  composite  taken  from  the  work  of  Tatarenko  et 
al15.  Measurements  of  insertion  loss  as  a  function  of  frequency  for  a  saturation  induction, 
47iM  of  1750  Gauss,  FMR  linewidth  of  1  Oe,  and  an  in  plane  applied  field  of  1700  Oe 
were  shown  to  vary  with  electric  field  by  over  100  MFIz  by  the  application  of  a  3kV/cm 
field  across  the  PZT.  A  schematic  of  the  device  as  input  to  FIFSS  is  shown  in  Figure  5. 

HFSS  Bandpass  Filter  /  Solid  Model 


z 


Figure  5.  HFSS  bandpass  filter  solid  model  from  Tarantenko  et  al. 

The  device  is  simulated  in  HFSS  as  a  symmetric  2-port  device  having  a  plane  wave 
input  to  port  1  which  excites  an  input  microstrip  line  (18x1  mm).  Radiation  couples 
through  the  GGG  substrate  and  into  the  YIG  region  (5mm  x  1.5  mm  x  110  pm).  It 
couples  under  FMR  back  through  the  GGG  substrate  and  out  into  the  second  microstrip 
line.  Microstrip  separation  is  2.5  mm.  Filter  transmission  (S21)  is  measured  as  the 
output  at  waveport  2  at  the  end  of  the  second  microstrip  line. 

Since  we  could  not  accommodate  the  magnetostriction  directly,  we  solve  static  cases 
for  varying  magnitude  internal  H  fields.  In  each  case  the  internal  H-field  is  uniform 
across  the  sample.  It  is  assumed  that  this  is  closest  equivalent  possible  to  simulating  the 
effect  of  the  ferrite  under  magnetostriction  at  this  time. 

Materials  properties  used  in  the  simulation  are  as  follows.  The  lower  substrate 
(alumina)  under  the  microstrip  lines  has  a  dielectric  constant  of  10  and  zero  loss.  To 
represent  a  lossy  dielectric,  a  dielectric  loss  tangent,  s"/s'  value  for  the  material  is 
applied.  The  greater  the  loss  tangent,  the  more  lossy  the  material.  The  GGG  substrate 
has  a  dielectric  constant  of  7.7  and  a  dielectric  loss  tangent  of  0.0002.  The  YIG  element 
has  a  dielectric  constant  of  12  and  a  dielectric  loss  tangent  of  0.0002,  and  a  magnetic 
saturation  value  of  1750  Gauss.  For  the  YIG,  AH,  the  full  resonance  line  width  at  half¬ 
maximum  relates  to  how  rapidly  a  precessional  mode  in  the  biased  ferrite  will  damp  out 


15  A.  S.  Tatarenko,  V.  Gheevarughese,  and  G.  Srinivasen,  Electronics  Lett.  42,  xx  (2006). 
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when  the  excitation  is  removed.  In  HFSS,  the  input  parameter  “delta  H”  varies  with 
frequency  and  is  therefore  calculated  within  a  simulation  by 

delta  H=AH  *  FREQ  /  Fmeas  (3) 

where  AFI  corresponds  to  the  measured  linewidth  at  frequency  Fmeas,  and  FREQ  is  the 
HFSS  swept  frequency.  In  the  simulations,  AH  =  3  Oe  at  10  GHz. 

The  magnetic  bias  direction  is  set  as  per  Tatarenko  et  al  to  be  parallel  to  the  sample 
and  perpendicular  to  the  microstrip  lines  and,  in  HFSS,  it  is  set  equal  in  magnitude  to  the 
magnetic  saturation  value  (1750  Gauss).  The  port  imepedance  for  this  device  is 
calculated  to  be  approximately  40  ohms. 

A  simulation  under  these  conditions,  with  H  internal  =  1000  Oe  results  in  the 
transmission  (S21),  reflection  (S 11)  and  power  loss  (dB)  as  shown  in  Figure  8.  The 
power  lost  is  computed  by  the  expression 

Plost(dB)  =  dB10[l-mag(Sll)2-mag(S21)2]  .  (4) 

Although  in  Figure  6  HFSS  appears  to  reproduce  the  ferromagnetic  resonance  peak  at 
9.37  GHz  expected  for  this  system  (with  a  loss  of  approximately  -3dB),  the  linewidth  is 
about  a  factor  of  10  greater  than  expected.  Since  the  linewidth  is  defined  as  3  Oe  at 
10  GHz  in  the  simulation  parameters,  we  would  expect  to  see  a  linewidth  of  this  order  in 
the  result. 

0.00 


-10.00 
S'  -20.00 

CD 

s 

£  -30.00 


-40.00 


-50.00 

8.5  9.0  9.5  10.0  10.5 

Frequency  (GHz) 

Figure  6.  HFSS  simulation  result  for  the  YIG  filter  device  schematically  represented  in  Figure  8. 
Active  YIG  sample  (5x1.5  mm  by  110  pm  thick).  FMR  resonance  is  observed  at  9.37  GHz.  The  loss 
is  approximately  -3  dB  at  9.37  GHz.  Hint  =  1000  Oe.  Magnetic  saturation  filed  =  1750  Gauss. 

In  order  to  now  attempt  to  model  the  effect  of  the  magnetostriction,  Figure  7  presents 
the  simulation  results  of  varying  the  internal  field  input  to  HFSS.  While  it  looks  like  at 
an  value  of  Hint=  2600  Oe  there  is  90  MHz  of  tuning  over  100  Oe,  it  was  unclear  whether 
this  represented  tunability  or  a  mode  hopping.  At  low  Hint  fields,  the  FMR  appears 
constant  at  around  9.3  GHz. 
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Figure  7.  HFSS  simulation  result  of  FMR  resonance  frequency  for  varying  internal  magnetic  field 
(Hint)  for  the  YIG  filter  device  schematically  represented  in  Figure  6.  Magnetic  saturation  field  = 
1750  Gauss.  AH  =  3  Oe  at  10  GHz. 

To  investigate  the  effect  of  mode  hopping,  simulations  were  performed  to  observe  the 
mode  structure  around  the  FMR  resonance  as  H;nt  passes  throughout  this  range. 
Simulations  commenced  at  Hint  =  2000  Oe  which  is  below  the  observed  resonance  region 
of  interest  in  Figure  7,  and  through  Hmt  =  2900,  at  which  point  the  negative  shift  in 
frequency  was  observed  to  be  a  maximum.  Finally,  we  calculate  at  Hint=3200  Oe,  where 
the  FMR  frequency  appears  to  return  towards  its  original  value.  Figures  8a)  -  e)  present 
the  results  of  these  simulations.  In  Figure  8(a)  and  8(b)  one  observes  the  primary  FMR 
transmission  resonance  shift  upward  in  frequency  from  9.37  to  9.47  GHz  with  increased 
Hint  as  expected.  Notice  also  in  each  case  the  presence  of  a  quasi-static  mode  at  about 
10.1  GHz.  In  the  next  two  Figures  8(c)  and  8(d)  we  observe  more  mode  mixing.  As  the 
internal  field  increases  to  2600  Oe,  the  transmission  resonance  shifts  to  9.57  GHz  and 
weakens  as  additional  modes  appear  at  lower  frequencies.  For  H;nt  increased  to  2900  Oe, 
a  small  resonance  remains  at  9.81  GHz  but  a  significant  amount  of  mode  mixing  occurs, 
particularly  around  the  quasi-static  resonance.  The  9.81  GHz  mode  becomes  very  weak 
and  noisy  as  we  observe  a  single  stronger  mode  entering  at  lower  frequency.  At  Hint  = 
3200  Oe  the  transmission  spectrum  returns  to  the  values  observed  at  lower  Hm[  fields.  It 
is  again  well  defined  and  the  FMR  appears  at  9.3  GHz.  Figures  8(c)  and  8(d)  illustrate 
that  there  is  significant  modal  structure  for  this  device  close  to  the  FMR  frequency  for 
fields  of  Him  in  the  range  2600  -  2900  Oe,  and  that  the  potential  for  mode  hopping  clearly 
exists.  It  would  be  necessary  to  optimize  the  device  design  to  reduce  the  number  of 
modes  present  in  the  FMR  tuning  region  in  order  to  observe  the  effective  tuning  more 
clearly. 
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Frequency  (GHz) 

Figure  8(a).  HFSS  simulation  result  for  the  YIG 
filter  device  schematically  represented  in  Figure  7. 
Active  YIG  sample  (5xxl.5  mm  by  110  pm  thick). 
FMR  resonance  is  observed  at  9.37  GHz.  The  loss 
is  approximately  -3.5dB  at  9.37  GHz.  Hint  =  2000 
Oe.  Magnetic  saturation  field  =  1750  Gauss. 


Frequency  (GHz) 

Figure  8(b).  HFSS  simulation  result  for  the  YIG 
filter  device  schematically  represented  in  Figure  7. 
Active  YIG  sample  (5x1.5  mm  by  110  pm  thick). 
FMR  resonance  is  observed  at  9.47  GHz.  The  loss 
is  approximately  -3.4dB  at  9.47  GHz.  Hint  =  2300 
Oe.  Magnetic  saturation  field  =  1750  Gauss. 


Frequency  (GHz) 

Figure  8(c).  HFSS  simulation  result  for  the  YIG 
filter  device  schematically  represented  in  Figure  7. 
Active  YIG  sample  (5x1.5  mm  by  110  pm  thick). 
FMR  resonance  is  observed  at  9.57  GHz.  The  loss 
is  approximately  -3.2dB  at  9.57  GHz.  Hint  =  2600 
Oe.  Magnetic  saturation  field  =  1750  Gauss.  Note 
the  second  resonance  at  10.11  GHz  and  -3.3  dB 
loss. 


Frequency  (GHz) 

Figure  8(d).  HFSS  simulation  result  for  the  YIG 
filter  device  schematically  represented  in  Figure  7. 
Active  YIG  sample  (5x1.5  mm  by  110  pm  thick). 
FMR  resonance  is  observed  at  9.26  GHz.  The  loss 
is  approximately  -3.45dB  at  9.26  GHz.  Hint  = 
2900  Oe.  Magnetic  saturation  field  =  1750  Gauss. 
Note  the  second  resonance  at  10.16  GHz  and  -3.5 
dB  loss.  Also,  a  weak  resonance  exists  at  9.81  GHz 
with  -4.2  dB  loss. 
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Figure  8(e).  HFSS  simulation  result  for  the  YIG  filter  device 
schematically  represented  in  Figure  7.  Active  YIG  sample  (5x1.5 
mm  by  110  pm  thick).  FMR  resonance  is  observed  at  9.27  GHz. 

The  loss  is  approximately  -3.dB  at  9.27  GHz.  Hint  =  3200  Oe. 

Magnetic  saturation  field  =  1750  Gauss.  Note  the  second  resonance 
back  at  10  GHz  and  -3.3  dB  loss.  The  system  has  returned  to  its 
original  state  as  in  Figure  8a). 

5.  TASK  3:  BENCHMARKING  PROTOTYPE  ME  DEVICES  AGAINST  SOA 

Data  compiled  in  Table  3  represents  the  input  solicited  from  engineers  working  in  the 
areas  of  tactical  radio  and  RADAR  when  queried  about  the  figures  of  merit,  present  and 
future,  for  key  system  components. 


Table  3.  Specifications  for  Components  used  for  Radio  and  Radar. 


Device 

Current  Materials 

Desired  Specs. 

Filters 

Piezoelectric  (SAW):  very 
small; 

Metal/dielectric:  (Stripline) 
simple  and  inexpensive  Qm 
=100-500; 

Quartz::  very  high  Qm  =  104 

20%  tunability 
around  center  freq; 

Qm  of  10,000  across 
bandwidth; 

<  -0.1  dB  loss; 

Isolators, 

circulators 

Ferrites:  handle  large  powers 
but  bulky  and  requires  large 
current 

Isolation  -lOOdB 
Control  complex 
impedance  <0.1  dB 
across  antenna 
dimensions  10/pc 

Phase 

shifters 

GaAs:  small  and  relatively 
inexpensive; 

Ferrites:  handle  large  powers 
but  bulky  and  requires  large 
current 

-0.1- -0.2  dB  RMS 
loss  across  all  phase 
and  freq.  settings 
Independent  control  of 
phase  and  amplitude 
lps  -  10ns  response 
time  $10  /  pc. 
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With  reference  to  Table  3,  and  the  assessment  perfonned  on  this  seedling  we  make  the 
following  conclusions: 

1 .  Single  phase  materials  still  need  significant  investment  to  understand  material 
physics  and  improve  material  properties. 

2.  Composite  materials  are  ready  for  straight  forward  implementation.  Sophisticated 
applications  require  additional  modeling  tools  and  materials  advances. 

3.  Composite  materials  are  better  suited  to  low  frequency  applications.  Millimeter 
wave  applications  have  needs  that  are  unmet  and  additional  investment  is  needed 
to  determine  if  multiferroic  based  composites  can  address  those  needs. 

4.  Multiferroics  materials  need  to  prove  their  superior  perfonnance  in  a  field 
crowded  with  many  technology  alternatives. 

6.  TASK  4:  RISK/BENEFIT  ANALYSIS  FOR  ME  DEVICES  IN  MILITARY 
SYSTEMS 

It  is  difficult  to  make  this  analysis  in  general  as  every  application  has  slightly 
different  needs.  Clearly,  the  field  of  ME  materials  and  devices  is  immature.  However, 
assuming  that  a  strategic  investment  will  result  in  successful  maturation,  the  easiest 
benefits  to  visualize  are  in  space-based  applications  where  reduction  in  component  size, 
weight  or  count  represents  a  significant  cost  savings.  For  example,  consider  the  non¬ 
linearity  of  a  traveling  wave  tube  power  amplifier  which  is  still  the  amplifier  of  choice 
for  high  power  satellite  communications.  Complicated  schemes  have  been  devised  to 
correct  this  nonlinearity  in  both  hardware  and  software  and  have  been  successful  for 
frequencies  up  to  K  band.  However,  as  both  military  and  commercial  applications  move 
to  ever  increasing  bandwidth  at  higher  transmit/receive  frequencies  new  strategies  must 
be  developed.  Table  4  gives  the  SWAP  metric  for  the  replacement  of  a  predistortion 
phase  equalizer  with  a  ferrite/piezoelectric  composite  which  is  directly  integrated  into  a 
waveguide.  The  integrated  composite  film-on-waveguide  approach  would  potentially 
substantially  improve  the  size  and  weight  over  currently  used  discrete  components. 
Though  the  power  requirement  is  marginally  higher,  the  cost  reduction  for  the  reduced 
weight  greater  than  x20. 
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Table  4.  SWAP  Metric  for  Film-on- Waveguide  Equalization. 


Solution  vs. 

Film- 

■on- 

W-band 

current 

waveguide 

(WR-10) 

waveguide 

SWAP  metric 

Individual 

device 

Platform 

Total 

Individual 

device 

Platform 

Total 

Size  (in3) 

2x1  O'5 

0.005 

0.044 

10.5 

Weight  (lb) 

0.0019 

0.5 

0.03 

7.5 

Power  (W) 

0.05 

12 

0 

0 

Device  cost  ($) 

~50 

~12K 

~1  K 

240K 

Launch  cost  ($) 

19 

5K 

312 

75K 

7.  SUMMARY 

The  interest  in  incorporating  ME  materials  in  electronic  devices  has  been 
reinvigorated  by  the  promising  electric-field-based  control  of  magnetization  or 
magnetization-based  control  of  polarization  in  monolithic  materials.  Primary  screening  of 
ME  monolithic  materials’  based  on  their  ferroelectric  and  magnetic  ordering  temperatures 
indicate  they  are  not  ready  for  application  in  devices  and  require  further  investment.  Thus 
attention  turned  toward  composite-based  materials  that  magnetoelectrically  couple 
through  a  strain  field,  which  have  known  and  beneficial  properties,  operate  at  room 
temperature  or  higher  and  can  be  modeled  in  a  quasi-static  fashion.  While  these 
composite  materials  may  exhibit  acceptable  tunability  and  loss  at  low  frequencies,  strain 
nonuniformities,  fabrication  sensitivity,  and  proper  evaluation  of  high  frequency 
performance  are  still  limitations.  HFSS  modeling  to  evaluate  high  frequency  perfonnance 
proved  difficult  because  it  cannot  account  for  non-uniform  magnetic  fields  nor 
magneto strictive  strain  effects.  Because  of  these  current  limitations,  we  recommend 
continued  funding  for  ME  materials’  research  and  for  development  of  modeling  and 
engineering  expertise  to  realize  the  potential  of  these  novel  materials. 
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